Direct experimental determination of the spontaneous polarization of GaN 
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We present a universal approach for determining the spontaneous polarization P sp of a wurtzite semicon- 
ductor from the emission energies of excitons bound to the different types of stacking faults in these crystals. 
Employing microphotoluminescence and cathodoluminescence spectroscopy, we observe emission lines from 
the intrinsic and extrinsic stacking faults in strain-free GaN microcrystals. By treating the polarization sheet 
charges associated with these stacking faults as a plate capacitor, P sp can be obtained from the observed tran- 
sition energies with no additional assumptions. Self-consistent Poisson-Schrodinger calculations, aided by the 
microscopic electrostatic potential computed using density-functional theory, lead to nearly identical values for 
P sp . Our recommended value for P sp of GaN is -0.022 ± 0.007 C/m 2 . 

PACS numbers: 73.63.Hs,73.22.-f,78.55.Cr,78.60.Hk,61.72.Nn 



Crystals with a singular polar axis are pyroelectric, i. e., 
they exhibit a spontaneous polarization P sp in equilibrium. 
An important class of materials with this property are semi- 
conductors that crystallize in the wurtzite (WZ) structure. 
Among these, GaN stands out as the material used for solid- 
state lighting 1 and power electronics. 2 While in the former 
case P sp affects both the color and the luminous efficacy of 
the light emitters, 3 it is exploited for transistor design in the 
latter example^ 

The physical origin of a non- vanishing P sp lies in the sin- 
gular polar axis of the WZ structure, while its magnitude is 
determined by the deviation from the ideal tetrahedral coor- 
dination of the atoms. In consequence, P sp critically depends 
on the ionic bonding contribution and the crystal field, in 
particular the difference between the internal cell-parameter 
u and its ideal value | 4 However, P sp is inaccessible for an 
infinite bulk crystal. The microscopic treatment of the polar- 
ization relies on the transformation between a reference state 
and the investigated system, 5 which means that only polar- 
ization differences are accessible, experimentally as well as 
theoretically. For P sp of WZ materials, the natural choice 
for such a reference is the zinc-blende (ZB) phase. 6,7 This 
phase has the same next neighbor configuration and a sim- 
ilar bond-length as WZ; 6 however, with P sp vanishing for 
symmetry reasons. 

A theoretical determination of P sp following this micro- 
scopic treatment requires methods of high sophistication 
such as density-functional theory (DFT). Two conceptual 
frameworks are used to define P sp as the difference from 
the ZB phase, namely, the electrostatic approach 8 and the 
Berry phase method. 5 Both DFT methods solve the refer- 
ence problem by considering WZ/ZB heterostructures to de- 
termine P sp . This approach has been implemented for vari- 
ous WZ semiconductors: BeO, 6 SiC, 8 and ZnO, 9 as well as 
the group Ill-nitrides 4,7 However, the accuracy of DFT is re- 
stricted as the value of u predicted by this technique depends 
on the choice of the exchange-correlation functional, with 
published values for GaN ranging from 0.3755 to 0.38154^ 
Unfortunately, P sp is far more difficult to obtain directly 



from experiment. The analysis of WZ heterostructures 
results in the difference of total polarization (including 
the piezoelectric polarization P pz ) between the constituent 
materials. 10 There has been only one attempt, based on an 
indirect thermodynamic approach, to arrive at an indepen- 
dent estimate of P S p- U Lacking a reliable experimental deter- 
mination of Psp, the values recommended in the literature 12 
are those obtained theoretically despite their acknowledged 
uncertainty. 

In this Rapid Communication, we deduce P sp of GaN from 
the emission energies of excitons bound to basal plane stack- 
ing faults (SFs) of the intrinsic Ii and I2 as well as the ex- 
trinsic E type. 13 These SFs are a local deviation from the 
hexagonal WZ (0001) stacking sequence to the cubic ZB 
(111) stacking sequence 14 with a structurally well-defined 
thickness. In other words, our samples are an experimental 
implementation of the WZ/ZB heterostructures used to the- 
oretically determine P sp . 6 Due to the smaller bandgap of the 
ZB modification, these SFs form perfect ZB-like quantum 
wells (QWs) in a WZ matrix, 15 ^ in the sense that fluctua- 
tions in composition and thickness do not exist. Moreover, 
the in-plane lattice constants of the two modifications are 
close to each other (Aa/a <2xl0 -3 ), 17 resulting in a neg- 
ligibly small contribution of P pz . In consequence, the inter- 
nal electrostatic field in the QWs formed by the SFs is di- 
rectly given by P sp in the WZ matrix. We use the observed 
differences in transition energies between the SF-types in 
conjunction with a parameter- free plate capacitor model to 
determine the strength of the field and thus P sp . These re- 
sults are quantitatively confirmed by self-consistent Poisson- 
Schrodinger calculations for different band alignments. The 
effective electronic width of the SFs needed for these cal- 
culations is obtained from the microscopic electrostatic po- 
tential computed by DFT within the superlattice approach. 
In this work, we realize the theoretical physical concept of 
WZ/ZB heterostructures 4-9 experimentally in order to quan- 
tify P sp for wurtzite semiconductors; for simplicity, we will 
refer to the value obtained following this definition as P sp of 
GaN. 
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FIG. 1. (Color) (a) /iPL spectrum of an individual microcrystal 
at low-excitation density. (D°,X) and FX are at the position of 
strain-free material, while individual SFs and SF bundles emerge as 
sharp lines right below the energies of the three SF types. The inset 
shows a top-view scanning electron micrograph of such microcrys- 
tals. (b), (c) False-color monochromatic CL maps superimposed 
on the corresponding cross-sectional scanning electron micrograph 
of a GaN microcrystal with the detected energies denoted in the 
respective colors. The CL related to SFs shows a characteristic 
elongation along the basal plane. 



The sample used in this study was obtained by pendeoepi- 
taxial overgrowth of self-induced GaN nanowires 18 using 
plasma-assisted molecular-beam epitaxy. The resulting GaN 
microcrystals are unique with respect to two important prop- 
erties: they are free of strain and dislocations and con- 
tain all three types of SFs. For the spectroscopic analysis 
of individual microcrystals, we employed a Gatan Mono- 
CL3 cathodoluminescence (CL) system and a Jobin-Yvon 
microphotoluminescence (/iPL) setup. The CL system is 
equipped with a photomultiplier and a charge-coupled de- 
vice (CCD) detector and mounted to a Zeiss Ultra55 field- 
emission scanning electron microscope. 19 To achieve a high 
spatial resolution (< 50 nm), the acceleration voltage was 
set to 2 or 3 kV, while the spectral resolution was chosen to 
be in the range of 1-5 meV. For /iPL, we used the 325 nm 
line of a Kimmon He-Cd laser for excitation 20 focused to a 
spot diameter of about 3 jim and attenuated to an intensity 
of 1 Wcm -2 . The spectral resolution was set to 1 meV. For 
all measurements, the samples were cooled to 10 K. 

Figure [Ha) presents a low-excitation /iPL spectrum of 
a single GaN microcrystal detached from the film and dis- 
persed onto a Si substrate. At high energies, emission from 




FIG. 2. (Color) High-resolution transmission electron micrographs 
revealing SFs of all three types and the bundling of such SFs with 
only thin WZ interlay ers. The micrographs (a) and (b) were taken 
on two crystals from the same sample used for the spectroscopic 
analysis. Both micrographs contain Ii SFs. In addition, micro- 
graph (a) contains the I2 and micrograph (b) the E SF. The overlaid 
ball- and- stick model (note that bright spots correspond to tunneling 
positions) illustrates the stacking sequence as also indicated explic- 
itly for each SF type. 

bound [(D°,X)] and free (FX) excitons is observed at the po- 
sitions expected for strain-free GaN. The spectrum is dom- 
inated by lower energy lines with a high-energy onset at 
about 3.42, 3.35 and 3.30 eV, which is associated with emis- 
sion from Ii , I2, and E SFs, respectively. Monochromatic CL 
imaging on the cross section of such a crystal, as depicted in 
Figs. [Hb) andQJc), shows that the near-band edge lumines- 
cence originates from the upper part of the microcrystals, 
while the remaining lines originate from the (lower) region 
of lateral expansion from a single nanowire to a microcrys- 
tal. All the low-energy lines appear as stripes running along 
the basal plane, confirming their common origin. 

Figure [2] shows high-resolution transmission electron mi- 
crographs of this transition region for two GaN microcrys- 
tals. All three types of SFs are observed in this region. The 
images also show that the SFs are often separated by only a 
few layers of WZ GaN. This "bundling" of SFs will modify 
the transition energy associated with a specific SF arrange- 
ment. For larger spatial separations, the electric fields within 
the aggregate of SFs will be redistributed such as to blueshift 
the transitions with respect to an isolated SF, 21 while a small 
separation will lead to a coupling of the associated electronic 
states, resulting in a redshift of the transitions. 22 Hence, the 
random separation of SFs of different types as observed in 
Fig. [2] leads to a (potentially) continuous distribution of the 
transition energies, explaining the multiple peaks observed 
in Fig.QJa). The transition energy of isolated SFs is, in con- 
trast, well-defined and is thus expected to dominate statisti- 
cally. 

Prior to compiling such a statistics, it is crucial to ensure 
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FIG. 3. (Color) CL spectral image acquired on the cross section of 
a microcrystal along the arrow indicated in the scanning electron 
micrograph in the inset. For each position along this line, a spec- 
trum is recorded which is displayed in the spectral image with its 
intensity color coded on a logarithmic scale. Note that the energy is 
blue- shifted by AE (26 meV for the marked case) due to a screen- 
ing of the polarization field when the electron beam approaches the 
SFs. 



that the internal electrostatic fields we are intending to mon- 
itor are not screened due to high excitation densities. The 
CL spectral image presented in Fig. [3] shows that screening 
may indeed occur: the crescent- shaped features reflect that 
the transition energy is lowest with the electron beam far off 
the SFs but steadily increases with the electron beam clos- 
ing in. The maximum blueshift is observed with the electron 
beam situated directly on the SF. The origin of this behavior 
is the increase in carrier density with increasing proximity 
of the electron beam due to diffusion and eventually direct 
excitation. 

To rule out excitation-induced shifts in the determination 
of the transition energies of the isolated SFs, we thus use 
exclusively (i) low-excitation /iPL spectra such as the one 
depicted in Fig. [Ha) and (ii) the low-energy onsets obtained 
for each of the transitions in CL spectral images. Figure 
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FIG. 4. (Color online) Histogram compiled of peak energies from 
low-excitation /iPL and of the lowest (unscreened) energies in CL 
(excluding phonon replicas). The Gaussian fit further illustrates the 
energies around which the peaks assigned to the different SF types 
are grouped. 



|4] shows a histogram of SF-related emission energies com- 
piled from spectra and spectral images recorded on several 
dozen microcrystals. The energies associated with the iso- 
lated SFs emerge in this statistical analysis of the peak en- 
ergies: we observe three peaks at 3.42 ±0.01, 3.35 ±0.01 
and 3.29 ± 0.015 eV in good agreement with the high-energy 
peaks for the three groups of lines in Fig.QJa). 23 In order of 
appearance, these peaks are associated with the Ii, I2, and E 
SFs, respectively, as expected from the thickness of the SFs 
(ZB segments) increasing in this order. While this confirms 
literature reports on the emission associated with the Ii and 
I2 SFs, 15 we also show that luminescence around 3.29 eV 
can be clearly attributed to excitons bound to E SFs. Emis- 
sion over the same energy range as displayed in Fig. |4] was 
recently reported for faulted GaN nanowires. 16 The obser- 
vation of emission below the band gap of bulk ZB GaN was 
attributed to the presence of internal electric fields, but the 
limited statistics did not allow for the identification of tran- 
sitions from isolated SFs and thus prevented an independent 
estimate of P sp . The accurate determination of the transition 
energies from all three types of SFs [cf. Fig.lH enables us to 
go further as detailed in the following. 

The crucial point in our analysis is that we rely on the 
differences in the emission energies rather than on their ab- 
solute values. The difference in thickness between each 
of the SFs is a single (111) bilayer (or molecular mono- 
layer) of ZB GaN, i. e., Ad = 0.259 nm. The differences 
in the emission energies are AEi^i 2 = 70 ± 15 meV and 
A£i 2 ^e = 60 ± 18 meV; thus the mean energy difference for 
adding a bilayer to the QW is AE = 65 ± 23 meV. Now let 
us assume that the spontaneous polarization is sufficiently 
strong so that the single particle energies (relative to the re- 
spective band) are governed by the triangular potential and 
thus remain basically the same for all SFs, while the change 
in confinement remains negligible. In this case, the polar- 
ization sheet charges at the interface between the ZB and 
WZ modifications essentially represent a plate capacitor, for 
which the addition of a slab of dielectric of the width Ad 
results in the potential difference AV = AE. From these ele- 
mentary considerations, we can directly calculate the polar- 
ization sheet charge density 

CT = |^ S p| = =0.021 ±0.007 C/m 2 , (1) 

where £ = 9.5 is the static dielectric constant for GaN and £0 
is the permittivity of free space. Note that the corresponding 
electric field within the SF amounts to 2.5 MV/cm. 

The same argument also applies to the microscopic elec- 
trostatic potential calculated in the present work and de- 
picted in Fig. 0a) for the E SF (for details, see Ref. 24 ). 
Tracing the triangular envelope of the potential in one super- 
cell [cf. Fig. [3a)], we estimate an energy difference AV = 
0.21 eV for a thickness of Jesf = 0-96 nm or 3.7 bilayers. 
Using Eq. this yields a value of P sp = -0.018 C/m 2 . The 
difference of this value compared to an earlier computation 4 
is due to the use of the projector- augmented wave method in 
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FIG. 5. (Color online) (a) Microscopic electrostatic potential V e s 
for the E SF calculated by DFT. From the envelope of the potential, 
we estimate the thickness (3.7 bilayers) and AV = 0.21 eV result- 
ing from the polarization field, (b) Band-profile for the E SF from 
a selfconsistent Poisson-Schrodinger calculation for a type-I band 
alignment. The horizontal lines denote the single particle energy 
levels, the arrow the resulting transition, (c) Transition energies 
for all three SF types from the Poisson-Schrodinger calculations 
with varied P sp and for two different conduction band offsets AEq 
corresponding to a type-I and a type-II alignment. The vertical 
lines denote the values of P sp for which the spectroscopically ob- 
served differences in the transition energy between the SF types are 
reproduced best; the absolute emission energies from experiment 
(including error bars) are shaded in gray for comparison. 



the present work, leading to a more accurate description of 
the electron density near the cores. 

To confirm the validity of the assumption inherent in 
the plate capacitor model used above, we performed self- 
consistent Poisson-Schrodinger calculations 25 for the QWs 
formed by the Ii, I2 and E SFs. We set the effective elec- 
tronic thickness of the E SF to 3.7 bilayers as motivated 
from the DFT results displayed in Fig. [3a). The thicknesses 
of the I2 and Ii SFs were then taken to be 2.7 and 1.7 bi- 
layers, respectively, in accordance with their differences in 
structural thickness. In order to examine the extent to which 
the transition energy depends on the respective band off- 
sets, we considered the two extremes suggested in the liter- 
ature, namely, AE C = 0.27 eV (Ref.QJ) and AE C = 0.15 eV 
(Ref.|24j), resulting in a type-II and a type-I band alignment, 
respectively. Standard values for the band gaps (containing 
a correction for the excitonic nature of the transition) 26 and 
effective masses were assumed. The residual donor density 
in our samples was set to Nd = 5 x 10 16 cm -3 . 27 Figure0b) 



TABLE I. Values for the spontaneous polarization P sp derived from 
theory and experiment. EM stands for effective mass and denotes 
the self-consistent Poisson-Schrodinger calculations. 



Method 



fip (C/m 2 ) 



DFT (Refs. 7, and i 24) 
Bond-orbital model (Ref. 29) 
Thermodynamics (Ref. 11) 
This work (DFT for E SF) 
This work (plate capacitor) 
This work (EM, type I) 
This work (EM, type II) 



-0.014 to -0.034 

-0.029 

-0.022 

-0.018 
-0.021 ±0.007 
-0.023 ±0.007 
-0.022 ±0.007 



shows an example of the band profile, the ground states, and 
the resulting transition energy for the E SF assuming a type-I 
band alignment^ 

Next, we varied P sp across the range of values found in 
the literature (see Table ID for both band alignments. The re- 
sulting transition energies for all three SF types are shown 
in Fig. He). Evidently, the limiting values either under- 
or overestimate the differences in transition energies. The 
observed AE = 65 ± 23 meV is best reproduced for P sp = 
—0.023 ± 0.007 C/m 2 in the case of a type I band-alignment 
and for P sp = -0.022 ± 0.007 C/m 2 in the case of a type-II 
band alignment. The fact that both band alignments result in 
essentially the same value for P sp signify that the assumption 
made in our initial analysis is valid: the influence of confine- 
ment is indeed small compared to that of the electric fields, 
and the SFs essentially behave as plate capacitors. 

All values for P sp obtained in this work and available in 
the literature are summarized in Table U Considering the 
uncertainty of the former, we recommend P sp = —0.022 ± 
0.007 C/m 2 . This value agrees exactly (probably fortu- 
itously) with the only other experimental estimation of P sp 
as reported in Ref. 1 1 and is also in fair agreement with that 
directly computed for the E SF in the present work (cf . Ta- 
ble©. 

We have demonstrated that the spectroscopic fingerprint 
of SFs in GaN allows the determination of the strength of 
the spontaneous polarization in an inherently parameter- free 
way. This approach may also be used to determine the spon- 
taneous polarization for other important wurtzite materials 
such as SiC and ZnO. 
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